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Abstract. In this paper wc survey some recent advances in the analysis of 
marginally outer trapped surfaces (MOTS). We begin with a systematic re- 
view of results by Schoen and Yau on Jang's equation and its relationship 
with MOTS. We then explain recent work on the existence, regularity, and 
properties of MOTS and discuss the consequences for the trapped region. We 
include an outlook with some directions for future research. 
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Given a Riemannian 3-manifold (M,g) and a symmetric (0,2)-tensor k on M, 
the triple (M, g, k) is an initial data set for flat Minkowski spacctime if and only if 
the overdetermined system of equations 

9ab = ffSf ~ D a uD b u 
(1) , ab D a D b u 



k ab 



(1 + D c uD c uy/ 2 

has a solution for some flat metric g . Here indices are raised using the metric 

9ab- 

This statement and its proof appear in the paper |31j by Pong-Soo Jang, who 
attributes it to Robert Geroch. The details of the calculation leading to ^ can be 
found in [3TJ Appendix], see also (TTJ Appendix C]. In his paper Jang sets out to 
generalize Geroch's approach to proving the positive mass theorem (based on the 
inverse mean curvature flow) from the case of time-symmetric initial data to the 
general case. 

Recall that the Geroch mass for a two-surface EcM with scalar curvature Rs 
and mean curvature H is defined by 



(2) 167T7« Gcroch (S) = y||L^(2R E -H|)d Ms . 

Date: June 23, 2010. 

' Research partially supported by Australian Research Council Discovery Grant #DP0987650 
and by the NSF grant DMS-0906038. 



1 



2 



L. ANDERSSON, M. EICHMAIR, AND J. METZGER 



(M,g,k) 



Figure 1. 



The Geroch mass is a specialization of the Hawking mass [35] to the time-symmetric 
case. The explicit form given here appeared in [35]. Although phrased slightly 
differently in [53], Geroch's argument for the positive mass theorem is based on 
the observation that mGeroch(S s ) is monotone increasing for a smooth family £ s 
moving in the normal direction with speed given by the inverse mean curvature, in a 
3-manifold of non-negative scalar curvature. This inverse mean curvature flow was 
later analyzed by Huisken and Illmanen [3D], who were able to prove monotonicity 
of the Geroch mass for a weak version of the flow, and use this to give a proof of 
the Riemannian Penrose inequality. 

In generalizing Geroch's argument, Jang introduces the defects 



(3) 



Slab(u) = g a b + D a uD b u 

t f , D a D b u 
k ab [u) = k ab - 



(1 + D c uD c ufl 2 



The condition that (M, g, k) forms initial data for Minkowski space is therefore 
equivalent to the condition that k(u) = and g(u) be flat for some function it. 
Taking the trace of k(u) with respect to the metric g(u) yields the quasilinear 
equation 



(4) 



D a uD 



1 + D c uD c u 



k a b 



D a D b u 



VI + D c uD c u 



= 0, 



which in particular must be satisfied for the height function of any spacelike asym- 
potically flat hypersurface in Minkowski space. This is Jang's equation. 

At this point it is convenient to note that g a b(u) is precisely the metric induced 
on the graph M = (x, u(x)) in the Riemannian product space (M xl,j + dt 2 ). Let 



~ab 

M ~ 3 - 



D a D b u 



1 + D c uD c uy/ 2 



tr^(fc) =g ab k ab . 

Then H A ^ is the mean curvature of the graph M, with respect to the downward 
pointing normal, and tr^ (fc) is the trace of the restriction to M of the pullback of 
k ab to the product M x K via the canonical projection ir : M x R — > M. Now we 
can write Jang's equation in the form 



(H) 



H^-tr^(fc) =0 
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(M.g.k) 

Figure 2. 

Assuming that the triple {M, g, k) is the induced geometric data for a hypersurface 
in a spacetime satisfying the dominant energy condition, the induced scalar cur- 
vature is non- negative modulo a divergence term (which of course can be large). 
Jang then, following the approach taken by Geroch in the case of non-negative 
scalar curvature, introduces a modified inverse mean curvature flow depending on 
a solution of Jang's equation, as well as an adapted Geroch mass that he shows to 
be formally monotone along his flow. If these steps outlined by Jang can be made 
rigorous, then his arguments lead to a proof of the positive energy theorem in this 
general situation. 

Jang's work has not been developed further due to the fact that an effective 
theory for existence and regularity of solutions of Jang's equation @ was lacking 
until the work of Schoen and Yau, who applied Jang's equation differently from the 
original intention by using it to reduce the space-time positive mass theorem to the 
time symmetric case. Further, it is not clear how to define an appropriate weak 
solution of the modified IMCF introduced by Jang. 

f .f. Jang's equation and positivity of mass. A complete proof of the positive 
mass theorem was first given by Schoen and Yau |5f j , for the special case of time- 
symmetric initial data. They then extended their result to general, asympotically 
flat initial data satisfying the dominant energy condition by using Jang's equation 
to "improve" the properties of the initial data in [55]. We describe here several 
aspects of Jang's equation which play a fundamental role in their work. 
Firstly, Jang's equation is closely analogous to the equation 

(5) H s +tr E (fc)=0 

defining marginally outer trapped surfaces E C M, where as above H E , trs(fc) are 
the mean curvature of E and the trace of k restricted to E, respectively. 

Equations of minimal surface type may have blow-up solutions on general do- 
mains and an important step in |55j is the analysis of the blow-up sets for the 
solutions of Jang's equation. At the boundaries of the blow-up sets, the graph of u 
is asymptotically vertical, asymptotic to cylinders over marginally outer (or inner) 
trapped surfaces - here the above mentioned relation of Jang's equation to the 
MOTS equation comes into play. 

Secondly, the induced geometry of the graph M of a solution of Jang's equation 
can be confomally changed to a metric with zero scalar curvature without increasing 
the mass. 

The fundamental reason for this is that the analogue of the stability operator for 
M , i.e., the linearization of Jang's equation, has, in a certain sense, non- negative 
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spectrum. Equation Q is translation invariant in the vertical direction. One of 
the consequences of this fact is that the non-negative lapse for the foliation of 
Mxl arising by this translation can be viewed as a principal eigenfunction of the 
linearization of equation (j4|), with eigenvalue zero. 

This spectral property allows one to prove the inequality 

(6) / R A - / 2 + 2|^ / 0| 2 >O 

Jm 

valid for any compactly supported Lipschitz function 4> on M, which in turn im- 
plies that the Yamabe invariant of M is non-negative. This means that (M, g) is 
conformal to a metric of non- negative scalar curvature. In fact, equation © is 
stronger, since the Yamabe operator has the factor 8 instead of the factor 2 in front 
of the \D^(f)\ 2 . The extra term in ((HJ) with respect to the Yamabe operator allows 
to control the change of the mass under this conformal deformation. We refer to 
section 14.11 for the actual calculation. In performing the conformal transformation 
to zero scalar curvature, the cylinders of the marginal boundary components are 
conformally blown down to (singular) points. 

The relationship between existence of solutions to Jang's equation, the existence 
of MOTS, and concentration of matter was observed and exploited in |56j . 

This problem has been revisited more recently by Yau [64] , and by Galloway and 
O'Murchadha [22]. 

In spite of a great deal of activity related to the positive mass theorem and 
minimal surfaces in the years following the Schoen-Yau and Witten proofs around 
1980, little attention has been paid to Jang's equation and MOTS from an analytical 
point of view until relatively recently. 

1.2. Existence and regularity of MOTS. As mentioned above, the close anal- 
ogy between Jang's equation, the MOTS equation, and the minimal surface equa- 
tion was exploited in the work of Schoen and Yau [55.. In particular, in that paper 
ideas from the regularity theory for minimal surfaces were applied to Jang's equa- 
tion. The positivity of the analogue of the stability operator, as discussed above, 
plays a central role here, completely analogous to the situation in minimal surface 
theory. The analogy between minimal surfaces and MOTS was further developed 
in [1] , where stability for MOTS was stated in terms of non- negativity of the (real 
part of the) spectrum of the stability operator. The positivity property of the sta- 
bility operator for a strictly stable MOTS was used there to prove local existence 
of apparent horizons. Further, in the curvature estimates for stable MOTS 
were developed along the same lines as the regularity estimates for stable minimal 
surfaces and for Jang's equation. 

MOTS are not known to be stationary for an elliptic variational problem on 
(M,g,k). This means that the direct method of the calculus of variation is not 
available to approach existence theory in parallel with minimal or constant mean 
curvature surfaces. The results in [53] lead Schoen [UJ to suggest to prove existence 
of MOTS between a trapped and an untrapped surface by forcing a blow-up of 
solutions of the (regularized) Jang's equation. In order to carry out this program 
one would like to construct a sequence of solutions to the Jang's equations whose 
boundary values diverge in the limit. The physically suggestive one-sided trapping 
assumptions proposed by Schoen are not sufficient to accomplish this directly. 

These technical difficulties were first overcome in [3] using a bending procedure 
for the data to convert the one-sided trapping assumption into the two-sided bound- 
ary curvature conditions necessary to solve the relevant Dirichlet problems, hence 
leading to a satisfying existence theory for closed MOTS, and subsequently in an 
independent approach using the Perron method in |17j . These two constructions 
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have established further features of MOTS related to their stability [U [3] , outward 
injectivity [3] in low dimensions, and almost-minimizing property |17j . These prop- 
erties confirm that MOTS are in many ways very similar to minimal surfaces and 
surfaces with prescribed mean curvature, which they generalize, even though they 
do not arise variationally except in special cases. 

The Perron method was used to solve the Plateau problem for MOTS in [T7] and 
also to extend the existence theory for closed MOTS to more general prescribed 
mean curvature surfaces that do not arise from a variational principle, including 
generalized apparent horizons (see in |18j . The combination of the almost min- 
imizing property of MOTS and the Schoen-Simon stability theory [35] introduced 
to this context in [T5] provide a convenient framework for the analysis of MOTS in 
arbitrary dimension, allowing techniques from geometric measure theory to enter 
despite the lack of a variational principle. 

1.3. Overview of this paper. In section[3]we introduce notation and give some 
technical preliminaries. Section |3J provides a systematic and detailed overview of 
the analysis of Jang's equation and the MOTS equation. As mentioned above, the 
solutions to Jang's equation in general exhibit blow-up, and boundaries of the blow- 
up regions are marginally outer (or inner) trapped. Section 13.51 explains how this 
fact can be exploited for proving the existence of MOTS in regions whose boundaries 
are trapped in an appropriate sense. Stability of MOTS is discussed in subsection 
13.61 where we also describe a new result on stability of solutions to the Plateau 
problem for the MOTS equation. Section [4] discusses in detail some of the main 
applications of Jang's equation in general relativity, including the positive mass 
theorem, formation of black holes due to condensation of matter and the existence 
of outermost MOTS. Finally, section [5] gives an overview of some open problems 
and potential new applications of Jang's equation and generalizations thereof. 

2. Preliminaries 

2.1. Initial data sets and MOTS. In this section we introduce the notation, 
sign conventions, and terminology used in this survey. Classical references for this 
material are [29] and |61j . 

An initial data set is a triple (M,g,k) where M is a complete 3-dimensional 
manifold, possibly with boundary, together with a positive definite metric g and a 
symmetric (0, 2) tensor k. In the context of general relativity, such triples arise as 
embedded spacclikc hypersurfaces of time-orientable Lorcntzian manifolds (M, g), 
referred to as the spacetimc, with induced metric g and (future directed) second 
fundamental form k. Hence, if 77 is a future directed normal vector field of M C M 
such that .9(77, 77) = -1, and if £, C £ T p M C T p M where p £ M, then k(£,C) = 
g{D^r]^Q. Here, D is the Lcvi-Civita connection of the spacetime. 

Now let E C M be an embedded two-sided 2-surface and let v be a unit normal 
vector field of £ C M. We write h for the second fundamental form of E with 
respect to v so that h(£, C) = g(D^u 7 Q for tangent vectors £, ( £ T p E for any 
p £ E, where D is the Levi-Civita connection of (M,g). We may think of v as 
a vector field along E C M so that I = v + n is a future directed null vector 
field of E when viewed as a spacelike 2-surface of the spacetime. Note that then 
x(£> := (h + C) = <j{D(l, C) for tangent vectors £, C £ T p T,. This symmetric 
(0, 2)-tensor \ is the null second fundamental form of E C M with respect to I. Its 
trace 6>s = trs \ is called the expansion of E with respect to the null- vector field I. 
Note that Qy, = He + trs k where Hs is the mean curvature of E C M with respect 
to the unit normal v. 

In many places in this survey, E C M has a clearly designated outward unit 
normal. If we compute the null-second fundamental form and expansion of such a 
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surface £ with respect to the corresponding future-directed outward null-normal, 
we will write Xy, ano - f° r emphasis. If £ is such a 2-surface, with 0^. = on E, 
then we say that £ is a marginally outer trapped surface or MOTS for short. 

If we use the future-directed inward unit normal to compute the expansion we 
write an d 0% ■ If 0£ = we say that E is a marginally inner trapped surface or 
MITS. For ease of exposition, we say that a two-sided surface £ C M is an apparent 
horizon if it is either a MOTS or a MITS, i.e., if H s + tr s (fc) = holds for one of 
the two possible consistent choices to compute the mean curvature scalar. 

The following spacetime analogue of the Bonnet-Myers theorem, whose proof 
uses the Raychaudhuri rather than the Riccati equation, lies at the heart of the 
Penrose-Hawking singularity theorems of general relativity (cf. Proposition 
4.4.3]): if the null energy condition holds in M, i.e. if for the spacetime Ricci 
tensor Ric(u, v) > holds for all points q £ M and null vectors v £ T q M, and if 
the expansion Oy, of £ C M with respect to the null vector field I along E is negative 
at p £ E, then a null geodesic emanating from p in direction l(p) has conjugate 
points within a finite affine distance. Physically, this means that the surface area 
of a shell of light emanating from E near p will go to zero before it has a chance to 
'escape to infinity'. 

Frequently, additional assumptions are included in the definition of an initial 
data set in the literature. In this survey we will always state such extra hypotheses 
explicitly when needed. Two such extra assumptions will be particularly relevant for 
us. First, recall [5TJ P- 219] that a spacetime (M,g) is said to satisfy the dominant 
energy condition if its stress-energy tensor T := Ric — 5R3 has the property that 
for every p £ M the vector dual to the one form — T(n, •) with respect to g is a 
future directed causal vector in T p M for every future directed causal vector r\ £ 
T p M. Here Ric and R respectively denote the spacetime Ricci tensor and spacetime 
scalar curvature. Note that this dominant energy condition implies the null energy 
condition used above in connection with the formation of caustics along light-like 
geodesies in the spacetime. If -q arises as above as the future directed normal vector 
field of a spacclikc hypcrsurfacc M C M, then one can use the Gauss and the 
Codazzi equations to express the normal-normal component (abbreviated by /j) 
and the normal-tangential part of T (written as a one form J) along M entirely in 
terms of the initial data (M, g, k). Explicitly, 

(7) i(R M -N 2 M + (tr M (fc)) 2 ) 

(8) div M (k - tx M {k)g) =: J 

where Rm is the scalar curvature of (M,g), D is its Levi-Civita connection, and 
where \k\ 2 M and trjv/(fc) arc the square length and trace of k with respect to g. 
Physically, /1 and J are, respectively, the energy density and the current density 
of an observer traveling with 4-velocity 77. A consequence of the dominant energy 
condition for the spacetime (M,g) is that \i > | J| must hold on M. By slight abuse 
of language one says that an initial data (M,g,p) satisfies the dominant energy 
condition if fi > \J\ holds. An important special case of the dominant energy 
condition is when M is a maximal slice of the spacetime, i.e. when trjvf(fe) = 0, 
so that the dominant energy condition implies that the scalar curvature Ra/ > 
is non-negative. In the case where k = 0, M C M is called a time symmetric slice 
or maybe more precisely a totally geodesic slice. In the case of time-symmetric 
data (M,g,k = 0), the dominant energy condition is equivalent to Rm > 0, and 
apparent horizons EcM are precisely minimal surfaces. 

For our discussion of the positive mass and positive energy theorem in subsection 
I4.1l we will also need the following definition (stated here as in [S5]): An initial data 
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set (M, g, k) is said to be asymptotically flat if there is a compact set K C M 
so that M \ K is diffeomorphic to a finite number of copies of K 3 \ B(0, 1) (each 
corresponding to an end), and such that under these diffcomorphisms 

\9ij ~ <%l + M\d P 9ij\ + \x\ 2 \d 2 9ij \ = Odx]- 1 ) and | R M | + \d p R M \ = M" 4 



J pqt 

as well as 



\hj\ + \x\\d p k t ,\ + \x\ 2 \d 2 pq k l3 \ = 0{\x\~ 2 ) and 



i=l 



0(|x|" 3 ) 



as := \J (x 1 ) 2 + (a; 2 ) 2 + (a; 3 ) 2 — > oo on each end. For the rigidity part of the 
positive mass theorem it will also be necessary to assume that the third and forth 
order derivatives of the metric are 0(|cc| -4 ). These conditions guarantee that the 
ADM energy 

Eadm = lib E l l=r (di9ij ~ dj9ii) jx\ dn2 

and the ADM linear momentum 

j=i \ x \ r 

are well-defined, see [H |51 S5]. The spacetime coordinate transformations which 
leave the asympototic conditions invariant are asymptotically Lorentz transfor- 
mations and the ADM 4-momentum vector -P^dm = (-^adm, ^adm) ^ s Lorentz 
covariant under such transformations. In particular, the ADM mass toadm = 
\f— P^ dm (Padm)p, is a coordinate independent quantity. Note that we will be 
sloppy in the sequel and refer to the positive energy theorem, i.e. the question 
whether -Eadm 0, by the term positive mass theorem, which would rather be 
appropriate for the statement wadm _: 0. 

2.2. Linearization of the expansion. Marginally outer trapped surfaces £ C 
M C M in general initial data sets are not known to arise as critical points (or 
indeed to occur as minimizers) of a standard variational problem described in terms 
of the data (M, g, k). (However, note that by definition a MOTS £ C M is a critical 
point for the area functional inside the future-directed null-cone of S C M.) In 
recent years, properties of MOTS akin to those of minimal surfaces have been 
introduced and investigated. In this subsection we discuss the linearization of the 
expansion 9^ with respect to variations in M. This sets the ground for a discussion 
of the natural notion of stability of MOTS in subsection 13.61 

Let E C M be a two-sided hypersurface and let ip s be a smooth family of dif- 
fcomorphisms of M parametrized by s G (— 6, S) so that ipo is the identity. Assume 
that 4- \ s =ofs = °~ + fv on E, where a £ r(TE) is a tangential vector field, where / 
is a smooth scalar function on E, and where v is a smooth unit normal vector field 
of E. The mean curvature scalar and tangential trace of k of E s := <^ S (E) can be 
viewed as functions on E via pullback by tp s , and they are smooth functions of s. 
Here we agree that the mean curvature scalar of E is computed as the tangential 
divergence divs(^) of v. It follows that 

(9) 

-(O e H s , ( 7> + --| s= o^H s , =-A E /-(|h|i J + Ric(i/ ) i/))/ ) and 



■ 

ds 

-(DvtTx(k),o-) + -^\ s=Q <p* s t^ m (k) = 2k(v,Dvf)+D v (ix M (k))f-(D v k)(v,u)f. 



ds 
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Here, As and Dy, denote respectively the non-positive Laplacian and the gradient 
operator of S with respect to the induced metric, \h\s is the length of the second 
fundamental form of E, and Ric and D are the ambient Ricci curvature tensor and 
covariant derivative operator. The literature often refers to the first identity as the 
Riccati or Jacobi equation, in particular when the tangential part a of the variation 
vanishes identically If X € T(TE) denotes the tangential part of the vector field 
dual to k(v,-) one can compute that 

{D u k){v, v) = -H s k{u, v) + (h, + (div M k)(v) - div s X. 

Using the Gauss equation and the expressions for the mass density /j in to 
substitute for Ric(V, v), and the current density J in (jHJ) to re-write divM k, it 
follows that 
(10) 

-(z? S E , ( T) + ^-| s=o ^:fe s 

as 

= -A E / + 2(X,D E /) 

+ (i R E -\\h + fc|| - J{y) — + div E X - \X\ 2 + ±0 E (2tr M (fc) - fe)) / 
=: W 

where 6*s s := Hs 3 +trs 3 (A:) is the expansion of E s . See [TJ (1)] and [231 (2-3)], and 
also [SSJ (2.25)] for an important special case of this formula that we will come 
back to in subsection 13.61 Following [23] we point out that when E is a MOTS and 
/ > 0, one can rearrange and express j Ls / more compactly as 

(11) div s (X - log/) -\X-Dv log/|| + i R E -\h + fc|| - J{v) - fi. 

Note that if the dominant energy condition fj, > \J\ is assumed, then the last three 
terms make a non-positive contribution to this expression. 

3. Analytical aspects of Jang's equation and MOTS 

3.1. Jang's equation. Given a compact subset fi C M with smooth boundary, 
F £ C 1 (57), and a function <fi defined on dfl we consider the quasi-linear elliptic 
partial differential equation expressed in local coordinates x = (x*,x ,x 3 ) as 

( H uV "\ ( D * D J U i \ v< \ n 
Q i ^ ,n — - hi = Fix) on fi 

(12) V l + \Du\y \Jl + \Du\2 



on dil 



where k = hjdx 1 (g) dx^ , g = gijdx 1 ® dx^ , giug = Sj, u l = g^djU are the compo- 
nents of the gradient Du = u l <9, of u, |Dm| 2 = g^diudjU is its length squared, and 
DiDjU = didjU — ^ijUk are the components of the second covariant derivative (the 
Hessian) of u so that Vdu = (DiDjUjdx 1 <£> (ix- 5 , where {Iy} are the Christoffcl 
symbols of the g metric in these coordinates. The left-hand side of (|12j) is indepen- 
dent of the choice of coordinate system and is easily seen to be of minimal surface 
type [551 Chapter 14]. Indeed, the graph of u viewed as a submanifold graph(u) of 
(M xl,j + dt 2 ) (where t is the vertical coordinate) and parametrized by an atlas 
of M via the defining map x — > (x,u(x)), has, in 'base coordinates', metric ten- 
sor gij + diiidjU with inverse g 1 ^ — , downward pointing unit normal vector 
v = 7 dj-dt anc j seconc i fundamental form ha = / P '' Pj " , so that the term 
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in (|12[) can be interpreted geometrically as the mean curvature of graph(u), while 
the remainder 

*(*)(«) := - Tq ^) % 

computes the trace of k (viewed as a tensor on the product manifold trivially 
extended, so that k(dt, •) = 0) over the tangent space of graph(w). We extend these 
geometric definitions in the obvious way to two-sided hypersurfaces £ C M X K 
and write Hs and trjj(fc) (they are functions on £). Note that in order to interpret 
the mean curvature scalar unambiguously we need to specify a normal vector field 
of S, while trs(fc) makes sense independently. If £ = graph(u) arises as above, 
Hs will always be computed with respect to the downward unit normal (i.e. as its 
tangential divergence) so that Hs(x, u(x)) = H(it)(x). When E is a hypersurface 
in the base M, then clearly trs(fc) = trsxR^c) and also He = He x r provided the 
orientations match. It is clear that, with appropriate identifications, the geometric 
operators H and tr(/c) are continuous with respect to C 2 and C 1 convergence of 
hypersurfaces respectively. 

3.2. Schoen and Yau's regularization of Jang's equation. In 53], R- Schoen 
and S.-T. Yau introduced the geometric perspective on solutions of Jang's equation 
discussed in subsection 13.11 and showed that solutions should only be expected to 
exist in a certain sense that we will explain in detail in subsection 13.41 They used 
their existence theory for Jang's equation to reduce the spacetime version of the 
positive energy theory (with general k satisfying the dominant energy condition) 
to the special case [51j of the positive energy theorem where the scalar curvature 
of the initial data set is non-negative. The analytic difficulty with Jang's equation 
H(u)— tr(fc)(u) = is the lack of an a priori estimate for sup n \u\ due to its zero order 
term tvM^k) = g lJ kij. That this is not just a technical difficulty but a fundamental 
aspect of Jang's equation will become apparent below in Theorems 14.51 14.31 The 
approach of |55| to bypass this issue is a positive capillarity regularization term: 

Theorem 3.1 (Schoen and Yau). Let (M,g,k) be an initial data set, f! C M a 
bounded subset with C 2,a boundary, <f> £ C 2,a (dfl), < r < 1, and assume that 
Hasi ~|trao(fc)| > T(f) along dft. Then there exists a unique solution u T G C 2 ' a (Cl) 
of 

H(u r ) — tr(fc)(u T ) = tu t on fi 
^ ^ u T = (f> on dQ,. 

If (M, <7, k) is asymptotically fiat ( cf. subsection \2.1\) . then there exists a unique 
solution u T e C 2 > a (M) of 

H(it T ) — tr(k)(u T ) = tu t on M with 

(14) 

u T —» as \x\ — > oo on each asymptotically flat end. 
Note that the estimates 

supt|u t | < max{3|fc| c( n),r|<?!>|c(9n)} 

respectively 

supr|u r | < 3|fc| c(M ) 

M 

are immediate from the maximum principle for solutions u T of these regularized 
equations (fT5|) . (|14[) . The second part of this theorem is proven in detail in [53] . 
The boundary gradient estimates necessary to establish existence in the first part 
can be derived from the boundary curvature condition Hgn — | trao(fc)| > Tip from 
a classical barrier construction due to J. Serrin (as described in [551 §14]): a suffi- 
ciently (C 2 -) small monotone inward perturbation of the boundary cylinder dQ, x M 
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below the rim {(x, 4>{x)) : x € dQ} will be the graph of a function u T defined near 
dfl which satisfies H(u T ) — tr(fc)(u r ) > tu t (because r > 0) with u T = (j> on 9f2 
and hence is a sub solution. The condition that Hgn + tr^n(fc) > t</> along dfl can 
be used to construct a super solution u T near <9f2 by perturbing inward above the 
boundary rim. For concise references and details see [T71 §2]. 

3.3. Classical results on minimal graphs and their limits. In order to moti- 
vate and explain the analysis [55] of the 'limits' of solutions as r \ of the regu- 
larized equations u T in Theorem 13. 11 we are going to review some classical results 
from minimal surface theory First, recall that if u <G C 2 (f2) satisfies the minimal 
surface equation, then E := graph(w) C O x R is an area minimizing boundary in 
Q x R. This means that if we denote by E = {(x,t) : x e f2 and t > u(x)} the 
super graph of u in VL x R, then 

V(E, W) < V(F, W) for every FcSlxl with EAF CC W CC !lxl. 

Here, V(F, W) denotes the perimeter of the set F in W (cf. P- 5]). This follows 
from a classical calibration argument using the closed 3- form 77 := (d vo\ g+ dt 2 ) \y 
where v = , (tfdi — dt) denotes the downward normal to £ thought of as a 

y/l+\Du\ 2 

unit normal vector field on 17 x R. Being an area minimizing hypersurface, E is a 
stable critical point of the area functional, and hence 

(15) / (|/i|! + Ric(^))0 2 < f \D s cj)\ 2 for all 

Je Je 

Here D-£(f> denotes the tangential gradient of </> along E, \h\s is the length of the 
second fundamental form, and Ric is the Ricci tensor of the ambient M x R. 
Using the Raylcigh quotient characterization of the first eigenvalue of an ellip- 
tic operator, one sees that the stability inequality p5[) is equivalent to the non- 
negativity of the Dirichlet spectrum of Ls on compact domains of E, where Ls / := 
— (As + Ric(i/, v) + / is the linearization of the mean curvature operator, cf. 
subsection 12.21 (with k = 0). We also remind the reader that stability is implied 
by the existence of a positive function / > on E such that Ljj / > 0. This can 
be seen either analytically, by use of the maximum principle, or more directly by 
integrating the pointwisc inequality 

< ^Ltf = (A s log / + |D E log /| 2 + |/i| 2 + Ric(i/, 1/)) 2 



over E. Then 



(16) 



2 + Ric(^, 1/)) 2 - f 2(Ds0,^slog/)-0 2 |Dslog/| 2 
JE 

follows from an integration by parts and an application of Young's inequality. The 
vector field —dt generates a family of diffeomorphisms tp s that act by downward 
translation. Clearly, the variations of E by Lp s also have zero mean curvature, 
and hence it follows from subsection 12.21 that Ls (\) = 0, where as usual v = 
yl + |-Dit| 2 - Hence ^ is a positive Jacobi field for Ls on E. Our point here is 
that the stability property of minimal graphs as expressed by (|15|) can be recovered 
without recourse to their stronger minimizing property We emphasize that this 
Jacobi field ^ measures the vertical component of the (upward) unit normal of E. 

Minimizing boundaries are smooth in ambient dimension < 7 and they form a 
closed subclass F of boundaries of locally finite perimeter in 51 x R with respect 
to current convergence (see e.g. [IB]). It follows that every sequence of minimal 
graphs has a smooth subsequential limit E.; — > E in F. If — denotes as above 
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the vertical part of the (upward) unit normal of £j, then Lj^ v i 1 = 0, and hence 
^S;^- < . holds where (i = |Ric|em) is a constant independent of i. This 
differential inequality has a non-parametric interpretation: the vertical part of the 
unit normal of the hypersurface £j is a non-negative super solution of a geometric 
homogeneous elliptic equation on £,; . This aspect of minimal graphs passes to their 
subsequential limit E. The Hopf maximum principle then implies that on every 
connected component of £ the vertical part of the unit normal either has a sign 
or vanishes identically. Put differently: subsequential limits of minimal graphs arc 
minimizing boundaries in f2 x K whose components are either graphical over open 
subsets of Q, or vertical cylinders whose cross-sections are minimizing boundaries 
in the base CI. This analysis of the limiting behavior of minimal graphs is carried 
out in detail in |41) . where concise references can be found. 

3.4. Behavior of graph(u T ) in the regularization limit. We consider a sequence 
of solutions {m t }o<t<i C C 2,a (f2) of H(it T ) — tr(fc)(u T ) = tu t as in Theorem 13. II 

For general k, their graphs £ T do not satisfy an apparent, useful variational 
principle. However, using the Jacobi equation to compute the variation of the 
mean curvature of £ T with respect to vertical translation, one obtains that 

(A St + Ric(l/ T , V T ) + \h r \l ) — =-v T (H(«r)) — 

where as before v T := yl + |£)u T | 2 , h T denotes the second fundamental form of 
£ T and |/i T |s T is its length, v T = -^-{u l T di — dt) is the downward unit normal of 
£ r , Ric is the Ricci tensor of (M xl,j| dt 2 ), and where we differentiate H(u T ) 
on the right as a function on the base f2. Expanding the right-hand side using the 
equation for u T , one obtains 

(A St + Ric(l/ T) V T ) + \h T \l) — < 10 (|fc|jlf IftrlE, + \Dk\ M ) —. 

Multiply this pointwise inequality by v T (f> 2 where <f> G C^(£ T ), integrate over £ T , 
and integrate by parts as in (fl6|) to see that there exists a constant j3 = P(\k\ c i(n)) 
so that 

f (m c {v T ,v T ) + \h T \l T )4> 2 < f \D^\ 2 + pj (\h T \^ + i)4> 2 

for all G C c 1 (£ r ). 

This almost looks like the stability inequality (|15p for minimal graphs. Schoen and 
Yau observed that this inequality, together with the a priori estimate for t\u t \ that 
follows from the maximum principle in the situation of Theorem 13. 11 the bound for 
H(u T ) hence resulting from the defining equation for u T , and a local area bound for 
£ T that follows from this bound on the mean curvature and a calibration argument 
as in subsection 13.31 can be used to derive pointwise curvature estimates for £ r as 
in jlH]- Moreover, they derived a geometric Harnack inequality for the vertical part 
^j- of the unit normal of £ T in the effective form sup Bpns ^ < 7infs pn s T ^- for 
any extrinsic ball B p centered on £ T and such that E>2 P C Q. X R. 

These results of |55j replace the variational arguments for minimal graphs in 
subsection 13.31 Hence, given a sequence {u r }o< T <i of solutions of (fTB")) so that 
supq T t\u t \ < oo, their graphs £ r have a smooth embedded subsequential limit £ 
as r \ 0. (The curvature estimates depend on the distance to the boundary of 
fl. Boundary barriers as in the proof of Theorem 13 . 1 1 can be used to show that £ T 
must remain bounded near dQ x R so that standard PDE techniques can be used 
to analyze the limit of u T near d£l.) From the Harnack estimate above it is then 
not difficult to sec that the components of the limit £ are complete, embedded, and 
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separated from one another by a positive distance, and that each is either graphical 
or cylindrical. 

Obviously, Hs — trs(fc) = must hold and it follows that the graphical compo- 
nents in the limit are solutions of Jang's equation, and that the cross-sections of 
the cylindrical components are apparent horizons of the base. These cross-sections 
inherit an orientation from the surfaces E r of the regularization limit, depending 
on whether these graphs 'blow-up' or 'blow-down' along them. The union of the 
graphical components of E is the graph of a function u defined on an open set 
Oo C 57, and this u must be unbounded on approach to <9S7o \ dd. Moreover, 
So smoothly asymptotes the boundary cylinder above the respective component 
of 90q \ 957, diverging either to positive or negative infinity, and it follows that 
dflo \ 557 is made up of apparent horizons. 

In summary, one has the following general existence result: 

Theorem 3.2 ([55)). Let (M,g,k) be a complete 3- dimensional initial data set, let 
57 C M be a bounded open subset with C 2 ' a boundary such that Hgn > |trgn(^)| 
holds along <957 where Hon is the mean curvature of 557 computed as the tangential 
divergence of the unit normal pointing out o/57, and where trao(fc) is the trace of 
k over the tangent space of <957, and let <f> € C 2,Q (<957). 

Then there exists an open subset Qq C 57 with embedded boundary consisting 
of <957 and of the union of two finite, possibly empty collections of smooth disjoint 
connected closed apparent horizons {E^~} and {E™}, as well as a function u <S 
C 2 < Q (<957U (l ) so that 

H(u) - tr(fc)(u) = on S7 
u = 4> on dfl 

u(x) —> +oo uniformly as dist(x, E^~) — > 0, and 

u(x) — > — oo uniformly as dist(x, EJ) — > 0. 

Computing the mean curvature scalar with respect to the unit normal that points 
into flo, the surfaces £T satisfy H s + +tr s +(fc) = while H s - — tr s -(fc) = 0. 

i i 3 3 

The same conclusion holds if (M, g, k) is asymptotically flat, Q = M , and where 
u(x) — > on each of the asymptotically fiat ends. 

3.5. Existence of MOTS due to blow-up. The analysis of [551 Proposition 4] 
outlined in the previous subsection contains the following, general result: given 
a sequence of functions {it T }o<r<i C C 2 ' a (Q) with H(it r ) — tr(k)(u T ) = tu t and 
sup n T t|u t | < oo, then for every O' CC fl there exists a sequence r, \ so that 
the graphs E Ti converge to a smooth hypersurfacc S C O' X 1. This limit S is an 
apparent horizon in the initial data set (M x R, g + dt 2 , k), and each component of E 
is either graphical or is a vertical cylinder whose cross-section is an apparent horizon 
in the base fl'. The sets f2± := {x € Q, : limsup^^ ±u Ti (x) = ±00} are disjoint 
and relatively open in 57', and their relative boundaries in 57' are smooth, embedded 
apparent horizons. Let 57o := 57\ (57 + U57_). The union of the graphical components 
of £ is a graph u : 57' n 57 — > R solving Jang's equation H(u) — tr(k)(u) = 0. The 
union of the cylindrical components of E is given by 57' n (<957 + n 957_) x R. The 
function u tends to ±00 near S7'n(<957on957±) and its graph is smoothly asymptotic 
as a submanifold to the vertical cylinder based on this set. In particular, if the 
sequence u T (x) diverges to +00 for some x G 57 as r \ while staying finite or 
diverging to —00 at other points of 57, then there must be (part of) an apparent 
horizon in 57. 

Using this blow-up analysis of Schoen and Yau and Theorem 13.11 for the ex- 
istence of solutions for the regularized Jang's equation, it is not difficult to see 
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that (non-empty!) closed apparent horizons E C fi exist in every bounded subset 
O C M, provided that dfl has at least two boundary components and satisfies 
Hgn — | tran(fc)| > (say greater than some small e > 0). Simply consider solutions 
u T : O — > R of H(u T ) — tv(k)(u T ) = tu t with u T = - on one of the boundary 
components and u T = — — on the others. From the maximum principle one has 
that sup^ t t\u t \ < max{e, 3|/c| C( -f2)}. The barriers used in the proof of Theorem 
13.11 can be used to show that these solutions u T diverge in a fixed neighborhood of 
the boundary. 

Schoen [37] suggested that the existence of MOTS should even follow if one only 
assumes that some boundary components of fi, whose union wc denote by 3+f2, 
satisfy Hd + n — trg + o(fc) > 0, while H^_o + tra_o(fc) > holds for the union cLSl of 
all other boundary components, because these conditions give rise to upper barriers 
for 'blow-up' near d+fl and to lower barriers for 'blow-down' near <9_f2. There is 
an important technical difficulty in implementing this approach: these one-sided 
barriers do not guarantee that solutions of the boundary value problems (|13p exist 
classically. 

That these technical difficulties can be overcome has been shown first by [3] and 
then by [17j using two independent methods which exhibit different features for the 
MOTS that are shown to exist. 

The 'bending of the boundary data' approach of [3] shows that the metric g 
and the second fundamental form tensor k of the initial data set {M, g, k) can be 
modified to g and A: in a neighborhood of d+Q. so that the region where the changes 
takes place is foliated by surfaces E s with positive expansion Hs„ — trs s (k) > but 
so that k vanishes identically near d+£l. A similar modification can be made near 
<9_£1. The above argument can then be used to show that a MOTS S exists in the 
modified data set. The maximum principle shows that this E cannot intersect the 
region where the data has been modified so that E is also a MOTS with respect 
to the original data set. The change of the geometry and the data are local but 
large. In particular, the change of the quantity |fc|c(n) that many of the geometric 
estimates for E in subsection !3.4l depend on is hard to control explicitly. Importantly 
though, [31 §4] develop a delicate barrier argument that shows that the MOTS E 
arising in the blow-up is stable in the sense of MOTS as discussed in subsection 
13.61 The curvature estimates of [2] , which do not require a priori area bounds and 
which depend only on the original data, are then available for E. 

An alternative line of proof is given in |17) . where the Perron method was 
introduced to the analysis of Jang's equation (cf. [55] for its classical applica- 
tion to minimal and constant mean curvature graphs) and used to generate max- 
imal interior solutions u T for the boundary value problems (|13[) with estimates 
supq jT t\u t \ < max{e, 3|fc|cm)}- These solutions won't assume particular bound- 
ary values in general, but they will lie above (respectively below) the lower (upper) 
barrier constructed from the boundary curvature conditions, which is all that is 
needed in the argument to force divergence near the boundary. Instead of using 
stability- based curvature estimates, a variant of the calibration argument mentioned 
in subsection 13.31 is applied to show that the surface E constructed in the process 
is a C-almost minimizing boundary in Q. By this we mean that E is the boundary 
of a set E in 51 such that 



(17) 



P{E, W) < V(F, W) + C\EAF\ 
for every Fc!l such that EAF CC W CC £1 



where C := 6\k\c(si)- (See [33] for a systematic study of such almost minimizing 
boundaries, and [TT] Appendix A] for further concise references on the relevant 
geometric measure theory.) So E minimizes area in £1 modulo a lower order bulk 
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term that is controlled explicitly. This feature of E C O is inherited from an 
analogous property of graph(u T ), see [TTJ Example A.l] for details. The stability- 
based curvature estimates in |55) are replaced by techniques from geometric measure 
theory which allow transitioning of the argument to high dimensional initial data 
sets if one accepts a thin singular set, as with minimal surfaces. 

In conclusion, we have the following existence theorem combining [5] and |17j : 

Theorem 3.3. Let (M,g, k) be a 3- dimensional initial data set and let C M be 
a connected bounded open subset with smooth embedded boundary dfl. Assume this 
boundary consists of two non-empty closed hypersurfaces d+fl and cLfi, possibly 
consisting of several components, so that 

(18) H-d+n - trg + o k > and Ha_Q + tr^_o k > 0, 

where the mean curvature scalar is computed as the tangential divergence of the 
unit normal vector field that is pointing out of fi. Then there exists a smooth 
closed embedded hypersurface Sell homologous to <9_f2 such that Hs +tr^(k) = 
(where Hs is computed with respect to the unit normal pointing towards d-Cl). E 
is stable in the sense of MOTS and it is C-almost minimizing in Q for an explicit 
constant C = C{\k\c<U))- This existence result and all the properties listed above 
carry over to initial data sets of dimensions < 7. In dimensions greater than 7 
we have the existence of a C-almost minimizing boundaries E in with a singular 
set of Hausdorff codimension at most 7 that satisfy the marginally outer trapped 
surface equation distributionally. 

The approach via the Perron method can easily be adapted to prove existence 
of surfaces E whose mean curvature is prescribed as a continuous function of po- 
sition and unit normal under boundary curvature conditions analogous to those 
in Theorem 13.31 see [TS], recovering classical existence results for variational pre- 
scribed mean curvature problems in special cases. In [17] , the Perron method has 
been used in conjunction with the analysis of Schoen and Yau described above to 
prove existence of MOTS spanning a given boundary curve, in analogy with the 
classical Plateau problem for minimal surfaces. We describe the general result for 
n-dimcnsional initial data sets: 

Theorem 3.4. Let (M n ,g,k) be an initial data set and let fi C M n be a bounded 
open domain with smooth boundary dQ. Let r n ~ 2 C dfl be a non-empty, smooth, 
closed, embedded submanifold that separates this boundary in the sense that dfl \ 
r"~ 2 = d-QUd+Q for disjoint, non-empty, and relatively open subsets <9_f2,9 + f2 
of dQ. Assume that Hgo — trg^ k > near <9 + f2 and that Hgn +tTon k > near 
<9_f2 where the mean curvature scalar is computed as the tangential divergence of 
the unit normal pointing out of f2. Then there exists an almost minimizing (rela- 
tive) boundary E™ -1 C f2, homologous to d-tt, with singular set strictly contained 
in 17 and of Hausdorff dimension < (n — 8), so that E Tl_1 satisfies the equation 
Hs+trj](A;) = distributionally, and so that E"" 1 is a smooth hypersurface near 
pn-2 boundary T n ~ 2 . In particular, if2<n<7, then E™ is a smooth 

embedded marginally outer trapped surface in ft which spans T n ~ 2 . 

We conclude this section by noting the close relation of the features of the regu- 
larization limit of Jang's equation from |55j with the classical Jenkins-Scrrin theory 
[33) of finding Scherk-type (i.e. infinite boundary value) minimal graphs in polygo- 
nal regions in R 2 and its obstructions, as expressed by the Jenkins- Serrin conditions, 
and also the extensions of this theory to infinite boundary value constant mean cur- 
vature graphs [S5] , [50] and [57] (see also references therein) in curvilinear domains 
in M 2 , S 2 and H 2 . 



JANG'S EQUATION AND MARGINAL SURFACES 



IS 



3.6. Stability of MOTS and an identity of Schoen and Yau. 

Definition 3.1 (|T],[2]). A closed two-sided surface E C M C M with vanishing 
expansion 9^ = ( computed with respect to the future- directed null normal I 
ii ■ u ■' r(S, TM) where v is a designated 'outward' unit normal vector field of 
E C M) is said to be stable in the sense of MOTS if there exists a positive function 
f > on E so that Ls / > 0, where 

L E / := -A S J + 2{X, D E /) + (I R s + fc|| - J(y) - fi + div s X - \X\ 2 ) f 

Here, X is the tangential part of the one form dual to k(v,-) on E and Re is the 
scalar curvature o/E. 

A few remarks are in order. First, note that L E / here is the linearization of the 
expansion 9 + for normal perturbations fv of E, cf. subsection 12.21 When k = 0, 
then this definition is consistent with the usual strong stability condition for closed 
minimal surfaces, as can be seen using the argument in subsection 13.31 Note that 
if we had / > with strict inequality L E / > 0, then a neighborhood of E in M 
could be foliated by surfaces {E s }_,5 <s<( 5 with Eo = E and so that E s lies 'outside' 
of E with respect to r\ and has positive expansion > for < s < 5, and such 
that E s lies inside of E and has negative expansion for — 6 < s < 0. (For minimal 
surfaces this strict stability condition implies that the surface is minimizing in this 
neighborhood.) Note that in general the operator Lg is not self-adjoint. It was 
noted in [TJ Lemma 1] that by the Krein-Rutman theorem the eigenvalue A of L E 
with the least real part is real, and that there exists an eigenfunction L E O = AO, 
positive on at least one and vanishing on all other connected components of E, 
corresponding to this principal eigenvalue A. The maximum principle then implies 
that the condition in Definition 13. II is equivalent to A > 0. As in subsection 12.21 it 
is useful to rewrite the pointwise condition < j L E / of Definition 13.11 as 

(19) M + J(u) < div s (A - D s log f)-\X- D s log /|| + i R s ~\h + 

If the dominant energy condition // > |J| holds, then the left-hand side here is 
non-negative, and an integration by parts exactly as in (|16p implies that 

(20) J l\ h + k\l< jT i R S + \D^\ 2 for every <f> G C^E). 
Together with the Gauss-equation one can conclude that 

(21) f \h\l<f< [ \D^\ 2 +P [ (|% + 1)0 2 

where j3 only depends on the initial data set but not on E. 

The following example of stable MOTS is one of the key observations in [5S] . Let 
u : floUdil graphical solution to Jang's equation in the sense of subsection 

13.41 So Slo C ri, the boundary of f2o consists of dil together with a finite number of 
smooth embedded apparent horizons, H(w) — tr(fc)(u) = on f^o, E := graph(u) C 
M x K is a complete submanifold with boundary {(x,u(x)) : x £ dil}, and u 
diverges uniformly on approach to the components of dilo \ dfl. As discussed in 
subsection l3.11 E has vanishing expansion in the initial data set (M x K, g + dt 2 , —k) 

(mind the sign!), and Ls / = where < f = —, — - — = = (—dt,u) is the normal 

yl+|i>u| 2 

component of the unit vector field generating downward translation. As above one 
has 

(22) n- J{v) =div s (X-D s log/)- \X- DsIog/|| + ^R^-\h-k\l. 

where /i and J are computed with respect to (M x R, g + dt 2 , k) (sorry!) and where 
X is the tangential part of the one form dual to —k(y, •) (we triple checked this 
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sign). Note that \i does not depend on the ^-coordinate and coincides with the mass 
density of (M,g, k), and that the same holds for J. This is equation (2.25) in [55] 
where it was derived by direct computation. See also equation (18) in |31| where 
the identity appears in disguised form and without geometric interpretation. 

It has been known from [SSJ that closed MOTS S C M that arise in the reg- 
ularization limit of Jang's equation are "symmetrized stable," i.e., the operator 
L^" 11 / := —As/ + (| Rs —\\h + fc|| — J(v) — fj) / on S has non-negative spec- 
trum. In [3] it was proven such surfaces are stable in the sense of MOTS, which is 
a stronger [211 Lemma 2.2] and physically more conclusive result. Here we discuss 
this stability from a geometric point of view, and we also discuss the stability of 
MOTS solving the Plateau problem in [T7j . 

Let E C M be a connected closed two-sided unstable MOTS with respect to the 
unit normal v. The Krein-Rutman theorem implies (cf. [T]) that there is A < 
and a strictly positive function O 6 C°°(E) so that Ls © = AO. The stability 
operator of ExR with respect to the extended initial data set {M xl,j + dt 2 , k) 

is LsxR = — jgz + Ls (where Ls ignores the dependence on the vertical variable t). 
Note that if T = T(t) is a function T G C 2 (R) then L SxK (OT) = 0(AT - J^T). 
Consider the function T(t) = — e (l — cxp(-^)) where e > is small and N > 1 is 
large. The relevant properties of T arc that T(0) = 0, that T < 0, that T{t) -e 
as t -> -oo, and that —T" + XT > > for all t G (-oo, 1] provided N 

is sufficiently large (depending only on A). Hence Ls x r(OT) > — O^ > t] in this 
range for a positive constant r\ > 0. This means that for s > sufficiently small, the 
surface {exp^ ^ (sO(9)T(t)v) : (8,t) 6Sx (— oo, 1)} is a smooth hypersurface (with 
boundary) in M x K that has positive expansion everywhere. Since T is monotone 
this hypersurface can be written as the graph of a function u : U — J- R where U is 
an open neighborhood of E such that u > in the part of U that lies to the side 
of v and so that u — > — oo on approach to the part of the boundary of U that lies 
in direction — v as seen from E. For u := —u we have that H(u) — tr(fc)(-u) < is a 
super solution of Jang's equation. Using — T instead of T one obtains a sub solution 
u of Jang's equation with analogous properties. (The awkward sign reversal here 
is due to the fact that Jang's equation is the MITS equation rather than a MOTS 
equation with respect to the data set (M xR,j + dt 2 , k).) 

There are three situations in which a closed MOTS E can arise in the regulariza- 
tion limit of Jang's equation: E C d£lo<~)d£l+, E C diloDdil-, and E C d£l-P\d£l+. 
The first two situations are the cases of graphical blow-up and graphical blow- down 
so there exists a solution of Jang's equation uq : Hq — > R which diverges to posi- 
tive, respectively negative infinity on approach to E C O^Iq. The strong maximum 
principle rules out the possibility that E be unstable in these cases straight away 
in view of the sub and super solutions constructed in the preceding paragraph. 
The third situation is the case of cylindrical blow-up: concretely, there is a family 
of graphs {ui}°l 1 C C°°(fl) where ut = u Ti solve the regularized Jang's equation 
H(it n ) — tr(/c)(u n ) = TiU Ti where \ 0, and such that it.; —X ±oo uniformly on 
compact subsets of fl±. From the analysis of Schoen and Yau in [SSJ Proposition 4] 
it follows that the hypcrsurfaccs graph(u;) c!!xM converge smoothly on compact 
sets to the 'marginally trapped cylinder' ExR. Note that graph(u^) is a super 
solution of Jang's equation where ui < and a sub solution where Ui > 0. Again 
using vertical translates of the sub and super solutions for Jang's equation above 
we can rule out the scenario that a component of a MOTS E arising in such a 
cylindrical blow-up be unstable. 

The next case to deal with is the Plateau problem. A MOTS E with boundary is 
stable in the sense of MOTS if there exists a function / on E, positive in the interior 
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and vanishing on the boundary, so that Ls / > 0. We note here that the Krein- 
Rutman theorem applies as before to show the existence of a real eigenfvmction 
(with Dirichlet boundary data) of Ls that is positive on at least one component of 
£ and vanishing on all the others. For MOTS with boundary we have the following: 

Lemma 3.1 f |19j). Assumptions as in Theorem \3.4\ in the smooth dimensions 
2 < n < 7. Then there exist solutions £ of the Plateau problem for MOTS in Q 
with boundary T that are stable in the sense of MOTS. 

Schocn and Yau used (|2"Tj) to derive pointwisc estimates for by adopting the 
iteration technique of |49j . They obtained the area bounds needed for this argument 
from a calibration argument for solutions of Jang's equation by comparison with 
extrinsic balls. These curvature estimates have been generalized in [5] to stable 
MOTS. This iteration technique extends to initial data sets (M,g, k) of dimension 
at most 6. It is remarkable and important that — as with minimal surfaces — in 
ambient dimension 3, stable MOTS have curvature estimates that are independent 
of a priori area bounds. This is uesed crucially in section f4. 31 

Theorem 3.5 (,2;). Let D, C M be a bounded open subset of an n-dimensional 
initial data set (M, g, k) where 3 < n < 6 and let E C O be a closed marginally 
outer trapped surface that is stable in the set of MOTS. Then one has the pointwisc 
bound 

(23) |/»| B <C(dist(E,5fi),|A!| CI(n) ,|mc w | c(n) ,iiij(n, ff ),|E|). 

When n = 3, then the bound on the right is independent of an a priori bound for 
the area | E | of E . 

We also mention that the regularity and compactness theory developed in [35] 
for stable minimal hypcrsurfaces generalizes to embedded MOTS, as was observed 
and used in |18j . This theory has the advantage of being available in all dimen- 
sions provided one accepts the usual singular set of Hausdorff co-dimension 7. This 
furnishes a convenient framework to carry out analysis on MOTS in high dimen- 
sions. This theory is particularly effective when combined with a one-sided almost 
minimizing property, see |18l Appendix A] and subsection 14.31 

4. Applications to general relativity 

In this section we discuss the applications that motivated the development of 
the mathematical theory for Jang's equation. 

4.1. The positive mass theorem. The first place where Jang's equation is ana- 
lyzed is in its application to reduce the general version of the positive mass theorem 
(PMT) to its time-symmetric form due to Schoen and Yau in [55] . 

The positive mass theorem is a question about asymptotically flat initial data 
sets (M, g, k) and its ADM-mass and ADM-momentum. 

Theorem 4.1 (Positive mass theorem). If(M,g,k) is a complete, asymptotically 
flat initial data set which satisfies the nn dominant energy condition fi > |J|, 
then mADM > |-Padm|- Moreover, if wiadm = 0, then (M,g,k) is initial data 
for Minkowski space. 

In the maximal case, where trj^f k — 0, the dominant energy condition implies 
R-m > 0. This leads to a formulation of the PMT relating only to the Riemannian 
manifold (M,g), called the Riemannian PMT. 

Theorem 4.2. Assume that (M,g) is asymptotically flat and has Km > 0. Then 
mADM > and equality holds if and only if (M, g) is flat K . 
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In a first step Schoen and Yau [5TJ[S3] proved the Riemannian PMT in dimension 
3 using the existence of certain area minimizing slices. Their method extends to 
dimensions 3 < n < 7 by a dimension reduction argument, see |52| and also |57j . 
The minimal surface argument of Schoen and Yau to prove the Riemannian PMT 
are closely related to their proof of the non-existence of metrics of positive scalar 
curvature on the torus in dimensions n < 7 in |50j . In fact, Lohkamp observed in 
[37] that the non-existence of such metrics and the time symmetric positive mass 
theorem are essentially equivalent in all dimensions. Two independent approaches 
to extend the positive mass theorem to all dimensions by addressing singularities 
of minimizing hypersurfaces when n > 7 have been given by Lohkamp |38| and by 
Schoen. 

An independent proof of Theorem l4 . 1 1 using spinor methods was later put forward 
by Witten [521 E2 • It does not need the reduction of the PMT to the Riemannian 
PMT that we describe below, and works in all dimensions under the topological 
assumption that the data set be spin. 

To describe the reduction of the general form of the positive mass theorem to 
the Riemannian case using Jang's equation, we follow [53]. The actual argument 
due to Schoen and Yau can be found in [53]. For the time being we assume that 
(M,g,k) is such that there exists a global solution u to Jang's equation (j4]) with 
boundary conditions u —> at infinity. By Theorem l3.2l we know that such solutions 
always exist provided M does not contain any closed apparent horizons. The graph 
M of u with the induced metric g is again asymptotically flat, and has the same 
ADM-mass as (M,g,k). The Schoen- Yau identity (|2"2"|) on M implies, in view of 
the dominant energy condition and a calculation similar to the one in section T3. 61 
that for all functions <j> £ C°° (M) with compact support 




Written in a slightly different way this implies 

(25) f 8|D^| 2 + 4? Rm > 6 f \D^\ 2 > 0. 

J M JM 

It then follows from standard methods that there exists a positive solution £ of the 
equation 

such that £ — > 1 at infinity. This implies that the conformal metric g := £ 4 g has 
scalar curvature R^ = 0. Moreover, it can be shown that £ has the asymptotic 
expansion 

C= l + A/r + 0(r- 2 ). 
Inserting £ as test function into equation (|25|) yields that 



A < - 




That £ is a legitimate test function can be verified by checking that the boundary 
term in the integration by parts, used to derive (f2~5j) from (f2"2"j) decays sufficiently 
fast. 

A direct calculation shows that the ADM-mass of (M, g) satisfies 

mxDu{M,g) = mxDu{M,g) + \A < m ADM (M,g) 

so that the resulting manifold (M, g) has ADM-mass no more than the initial data 
set (M, g, k). Since the scalar curvature is zero, the Riemannian PMT gives that 
mADM{M, g) > 0. If itiadm(M 7 g) = one can work backwards through this 
argument to see that in this case g is flat, £ is constant, so that g = g. Moreover, 
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there is also equality in the Schoen-Yau identity, so that h = k. In particular, the 
criterion of Jang ([1]) is satisfied and (M, g, k) is a data set for Minkowski space. 

Recall the simplifying assumption that a solution to Jang's equation exists on 
(M,g,k). This is indeed a restriction, as Jang's equation can blow-up (or down) 
asymptotic to cylinders over marginally outer (or inner) trapped surfaces, cf. The- 
orem 13.21 The resolution of the situation was achieved in |55j by compactifying 
the resulting cylindrical ends using a conformal transformation. While the actual 
procedure is out of the scope of this article, we wish to remark that this step is a 
major obstacle in the reduction of the general Penrose conjecture to the Ricmann- 
ian version, proved by Huisken and Ilmanen |30j and Bray [3] • A detailed discussion 
of this fact can be found in [39:. 

This reduction has been described by Schoen and Yau |55j for three dimensional 
initial data sets. The technical difficulties in higher dimensions are due to the 
potential singularities of apparent horizons and hence the blow-up cylinders in 
the solutions of Jang's equation, and their potentially complicated topology that 
prevents direct application of the arguments from |55j . These technical difficulties 
are resolved in |16j in dimensions 4 < n < 7. 

4.2. Formation of black holes. The mechanism that causes Jang's equation to 
possibly blow-up along apparent horizons yields an approach to the existence of 
apparent horizons in the following way. Assume that (M, g, k) is an initial data 
set, where M is compact with non-empty boundary. In addition suppose that the 
boundary geometry is such that the barriers needed to solve Jang's equation exist. 
Then the condition that (M,g,k) does not contain apparent horizons implies that 
the Dirichlct problem to Jang's equation is solvable without the possibility of blow- 
up with arbitrary boundary data, cf. Theorem 13.21 If one can devise conditions 
that lead to a contradiction using this solution, the existence of apparent horizons 
can be concluded. 

The first time that this prototype was used, is in the paper by Schoen and Yau 
to prove the following theorem. 

Theorem 4.3 ([5B]). Let (M,g,k) be a compact initial data set with non-empty 
boundary DM such that Hqm > | trgjvf (&)|- Let C M such that the following 
conditions are satisfied: 



Then M contains an apparent horizon. 

Here Rad(f2) denotes the H-radius of a set il which is defined as follows. Let 
r C SI be a curve bounding a disk in il. The radius of f2 relative to T is defined as 

Rad(fi,T) := sup{r : dist(r,<9Q) > r, T does not bound a disk in T r (T)}. 

Here T r (T) is the tubular neighborhood of V with radius r. The radius of then 
is defined as 



Roughly speaking, Rad(fi) is the diameter of the largest tubular neighborhood of 
a curve T that does not contain a disk spanned by T. 

As already indicated the argument proceeds via contradiction, so assume that 
there are no apparent horizons in M. Then the Dirichlet problem for Jang's equa- 
tion on (M, g, k) is solvable with zero boundary values, cf. Theorem l3.2l We denote 
the graph of the solution by M. From the Schoen-Yau identity (|22p , it follows that 




Rad(f2) = sup{Rad(i7, r) : T bounds a disk in ft}. 
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on the portion Mq of M above one has 

R A y > 2A + 2\lo\ 2 - 2 div^ us, 

where lu = X — D A[ log /, / = —(dt, 0), v is the downward unit normal to M, and 
X is the tangential part of the one form —k(v, ■) as before. This inequality yields 
that the first Dirichlet eigenvalue A of the operator L := — + | R A ^ on Mq 
satisfies A > A. Furthermore, since the distances in the g metric are no less than 

in the g metric, it also follows that Rad(Mn) > \f~^K 7T - 

The first Dirichlet cigenfunction of L on Mq, satisfying — A^<f)+ | <f) = \<fi 
is positive and can be used to define the following functional on surfaces E C M, 

A (E) = / cj)da, 

where da denotes the area element induced by g. Note, that this functional can 
be interpreted as the area functional for surfaces of the form E x S 1 in M x S 1 
equipped with the metric g = g + (f> 2 ds 2 , where ds 2 denotes the standard metric 
on S 1 . Note that g has scalar curvature R = R^ — 2tfi~ 1 A A j<f) = 2A > 2A in 
Mq. This interpretation also implies that one can find a minimizing disc E for 
with boundary T, where T is chosen such that Rad(Mn,r) > Rad(Mfj) — e. The 
minimizer E satisfies the Euler-Lagrange equation H = — {D^ log <f>, v) , where v is 
the normal vector field on E used to define H. More importantly, stability of E 
implies that the operator defined by 

/ h+ -A s / - (D E log <j>, Ds/) + /(A - \ R s +<F X &2<t>) 

has non-negative Dirichlet spectrum. The form of this operator follows for example 
by reduction of the stability operator in (M x S 1 ,^) for equivariant variations on 
surfaces with S' 1 -symmetry. Let ip > be the first eigcnfunction of this operator 
and define a functional for curves 7 in E as follows: 



-W(7) = 




Recall Bonnet's theorem, which asserts that stable geodesies in surfaces with scalar 
curvature bounded below by a positive constant have bounded length. Here, the 
modification of the length functional by introducing the weight (jjijj into 1^ leads 
to a similar effect for curves minimizing 1^ in the sense that the stability of the 
minimizer forces the minimizer to be short, in particular the length is bounded by 

w. By definition of Rad(M, T) the disc E intersects the boundary of a tubular 

neighborhood T r (T) of radius r < Rad(M, T). Thus it is always possible to find 
a minimizer for 1^ with length at least r, since one can minimize 1^ among all 
curves with one endpoint on T and one endpoint on dT r (T) n E. This yields the 
desired estimate for Rad(M), and thus for Rad(M). See [SS] for details. 

There are several variations on this theme in the literature. Clarke |13j gave 
an interesting and useful observer independent condition on the energy-momentum 
tensor of a space-time that implies the trapping condition on the boundary of the 
initial data set in Theorem 14.31 

A refined criterion for the existence of horizons was given by Yau |64j . 

Theorem 4.4. Let [M, g, k) be an initial data set satisfying the following condi- 
tions: 

(1) There exists c > such that Hqm — \ tr^M (k)\ > c. 

(2) Rad(M) > Jl^ir where A < |c 2 + /j - \ J\ on M. 
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Then M contains an apparent horizon. 

It is instructive to consider the case k = first: if c > is large, the first 
condition suggests that M shrinks quickly from its boundary dM inwards, while 
the second condition implies that the interior of M is large in a certain sense. The 
conclusion is that part of the interior of M must be separated from the boundary 
by a minimal surface. 

Again, Jang's equation enters prominently. Yau's argument in j64j is by con- 
tradiction and proceeds as follows. Assume in virtue of Theorem 13.21 that Jang's 
equation has a global solution u on M. Denote the graph of u in M x R by M and 
its induced scalar curvature by R M . Then, by the Schoen-Yau identity (|22p one 
has that 

2(// - | J|) < R Kl -2M 2 + 2div M u>, 
where oj = X + D M logv as in equation (f2"2"j) . where we use v to denote This 
yields for all <j) <E C°°(M) the following estimate: 

(26) 2 f O-|J|)0 2 < / 2\D M cf > \ 2 + R M cf > 2 + 2 f cj> 2 (uj,N), 

JM JM JdM 

where N denotes the outward pointing normal to dM in M. The point is that 
the difference of the boundary term in equation (|26[) and the mean curvature of 
the boundary has a positive lower bound, as one can see as follows. Recall that 
(lj,N) = (D M logv, N) - k(i>,N) where v = f^ 1 = sjl + \D M u\ 2 . Moreover, a 
calculation shows that the mean curvature of dM in M satisfies R dM = v^ 1 HdM, 
where the latter is calculated with respect to the metric g. The normal N is 
given by N = v^ 1 ^ + \D^ju\dt), where rj is the outward pointing normal to dM 
in M. To calculate Hg^ — {lo, N), note that since u — on dM we have that 
Dmu = o-\Dmu\t], where a G {±1}- We let V — 7r»z> = v _1 DmU = crv~ 1 \DMu\i r i. 
Then the mean curvature of M on dM is given by 

H A/ = divM V = adiv M (v~ 1 \D M u\n) = av^ 1 \D M u\ B dM +v~ 3 D 2 M u(n,r]) 

= ay- 1 1 D M u | Rom +<J \ D M u \ - 1 (D Kl log v, N). 

Note that this is also true if Dmu = 0, since then also Dm log u = 0. By Jang's 
equation, H A/ = tr fl - / (fc), where 

tr Af( fc ) = tr A/(fc) - k(i>, v) = tr dM (k) + v~ 2 k(r], rj). 
Since furthermore k(i>,N) = av~ 2 \D Mu\k{r),rj) it follows that 
= <r\D M v\(R & -tr A (k)) 
= v-^Dmu] 2 Rqm +{D m logv, N) - ati 9M {k) - k{v, N) 

and thus 

(u,N) =atx dM {k) - \D M u\ 2 R dM ■ 

Finally, 

^dM^^^) = " H 3M -<j\D m u\ tv dM {k) > v(R dM -\ tx d M{k)\) > c. 

This boundary term then has a similar effect as the A in an extension of the argu- 
ment of Schoen and Yau to get an estimate on Rad(M) contradicting the assump- 
tion as before. 

Of the several different proposals to define the size of a body in an alternative 
way, we want to mention specifically the suggestion of Galloway and O'Murchadha 
j2"2"] . They use the intrinsic diameter of the largest stable MOTS bounded by curves 
in the boundary of the body to define the radius of the body and show in turn that 
this radius is bounded if the matter content of the body is large. 



22 



L. ANDERSSON, M. EICHMAIR, AND J. METZGER 



The boundary effect discovered by Yau plays a crucial role in the proof that the 
quasi- local mass defined by Liu and Yau is non- negative [351 I36j . The common 
theme with section 14.11 is that Jang's equation is used to transform the question 
whether the Liu- Yau mass is non-negative to a question in Riemannian geometry. 
As before the transition to Jang's graph is followed by a conformal transformation 
to a metric with zero scalar curvature. In the Riemannian setting established by 
this procedure, the Liu- Yau mass is transformed to a quantity bounded below by a 
modified version of the Brown- York mass. This uses the boundary effect calculated 
above in a crucial way. Liu and Yau show that this quantity is non-negative by 
extending an argument of Shi and Tarn [55] , 

Eardley [TS] uses Jang's equation to give a different criterion for the formation of 
black holes. To this end, for a data set (M,g, k) and a region f2 C M the following 
quantity is introduced 



Note that /c m ; n is the smallest value that trfc(u) could take at any point of f2 for 
any graph u : — >• R. In the setting below, where k is positive definite, it equals 
the minimal value of the sum of two smallest eigenvalues of k on Q. 

Theorem 4.5 ([15]). Given compact initial data (M, g, k) with non-empty boundary 
such that > | i^dM(k)\- If there is O C M such that 

k min {Q)Vol(Q) > Area(dD), 

then there is an apparent horizon in M . 

Proof. The proof of this theorem is by contradiction. If there are no apparent 
horizons in M, then there exists a global solution of Jang's equation to the Dirichlet 
problem with zero boundary data. Denote the graph of this solution by M and by 
V = 7T* i>, the orthogonal projection of the downward unit normal. Then Jang's 
equation is equivalent to 



since db/M V is the mean curvature of M with respect to v and the right hand side 
is just the trace of k on M. Integrating this on 17 C M yields 



Here v denotes the outward normal to dQ in M . This contradicts the assumptions 



4.3. Existence and properties of outermost MOTS. In section |4~T1 the poten- 
tial blow-up of Jang's equation at apparent horizons is an undesirable property that 
has to be overcome. In section 14.21 the existence of apparent horizons is a rather 
indirect consequence. In contrast, the way Jang's equation is used to construct 
MOTS in section 1^751 is far more direct and can be used to derive crucial properties 
of outermost MOTS. 

To get started, fix a complete initial data set (M, g, k), and assume for simplicity 
that M is compact and that dM satisfies 0^ M > 0. We say that a MOTS E C M 
is outermost if it is of the form E = 90, where O C M, and the following holds: If 
E' = 90' is any other MOTS, with Q.' D fi, then 0' = £1 In other words, if E is 
outermost, then there is no MOTS in the region M\Cl exterior to E. 

We expect the outermost MOTS to be the boundary of the trapped region. To 
this end, we define a set C M to be trapped, if 9qq < 0. The trapped region T is 



k min (n) :=M{ti M (k)-k(v,v) :pen,veT p M,\v\ < 1}. 



dW M V = tr M (k)-k(V,V), 




of the theorem. 



□ 
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then the union of all trapped sets [HJ [53] , 

t= |J n. 

Q is trapped 

Using a slight extension of the existence Theorem 13.31 adapted to weakly trapped 
boundaries, cf. [3J Section 5] and also [TSl Remark 4.1], it follows that a trapped 
region f2 as in the definition of T is contained in a trapped region f2' D f2 whose 
boundary 90' is a MOTS, and such that d£l' is stable in the sense of MOTS and 
is C-almost minimizing with respect to variations in M . 

To conclude smoothness of dT as for example in [30j where the time-symmetric 
case k = is discussed, we need to verify three points. These are whether two 
intersecting MOTS are contained inside one smooth MOTS that encloses them, the 
embeddedness of dT, and area bounds. 

The question whether the union of two trapped sets is a trapped set relates to 
the following problem. Given a sequence of MOTS, E n = <90„, we wish to replace 
it by an increasing sequence T,' n = dfl' n so that fl' m C $l' n for all m < n, as in [50] . 
This can be handled in two different ways. In [3] a sewing lemma due to Kriclc and 
Hayward [34j was employed in conjunction with Theorem 13.31 to conclude that if 
two sets Hi and f^2 with 9$ n . = intersect, then there is f2 D OiUf^ with 9^ = 0. 
Alternatively, the Perron method and an approximation argument can be used to 
find an enclosing MOTS PH Remark 4.1]. 

To conclude embeddedness of dT we have to show that the limit of such an in- 
creasing sequence of MOTS E„ = dQ n is embedded. Since all the £„ are increasing 
and embedded, the only crucial point is that the limit £ does not touch itself on the 
outside. For minimal surfaces this scenario would be ruled out by the maximum 
principle, which does not work for MOTS in this situation. The problem is that 
locally two sheets of a MOTS may touch, but with opposite orientation. The case 
of two touching spheres in flat space illustrates this. To show that this can be ruled 
out for outermost MOTS, in [3] a quantity called the outward injectivity radius was 
introduced. It is then shown that one can assume it to be bounded below along 
the sequence £„ as above. This bound yields a lower bound on the arc length of a 
geodesic starting on £ n in direction of the outer normal, before it can intersect £„ 
again. The argument in [3] derives this property from the fact that along a short 
geodesic that joins two points on E n , a neck with negative 9 + can be inserted. 
Then the sewing lemma can be used to produce a barrier suitable for Theorem l3.3l 
This procedure can only be applied a finite number of times, since it can be shown 
that each surgery can be made at a place where it consumes a fixed amount of vol- 
ume outside of the initial MOTS. This surgery requires curvature bounds for stable 
MOTS, which have been derived in [2] in ambient dimension 3. Alternatively, one 
can use results from [IS] based on the lower order properties of horizons and the 
regularity theory of Schoen-Simon to conclude embeddedness of £. In fact, it is 
easy to see that if two sheets of the hypersurface £„ are close on the outside as 
above, then they can be joined by a small catenoidal neck to save area. This would 
contradict the almost minimizing property of £„ with respect to variations in the 
complement of fi n . This approach also works in higher dimensions. 

This leaves as a last point the fact that the area of the surfaces £„ needs to be 
bounded. In [18j these bounds are immediate from the almost minimizing property. 
In [3] it is shown that a lower bound on the outward injectivity radius implies an 
upper bound on the area. This follows from the observation that, given curvature 
estimates, the area of the MOTS can be estimated by the volume of the outward 
part of an embedded tubular neighborhood of radius p divided by p. Since the 
outward injectivity radius is bounded below for outermost MOTS, one can take a 
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fixed p and conclude the area bounds from the fact that there is only finite volume 
outside the MOTS. The approach in [3j is specific to ambient dimension three, since 
the lower bound on the outward injectivity radius requires the surface in question 
to have curvature bounded independently of the area. 

Let us investigate the topology of the outermost MOTS. In three dimensions 
outermost MOTS (assuming an outer untrapped barrier) are unions of topological 
spheres. This is well known in the time symmetric case in three dimensions, where 
MOTS are minimal surfaces, for example [20j uses minimal surfaces techniques from 
[42) or [301 Lemma 4.1] where this is proven without curvature restriction. 

For MOTS the question of topology was answered by Galloway and Schoen [33] , 
who showed that any stable MOTS must be of non-negative Yamabe type, provided 
the dominant energy condition holds. The argument is based on the Schoen- Yau 
identity, which follows for stable MOTS and a calculation similar to section 14.11 
Galloway [21] was able to exclude the marginal case for smooth outermost MOTS. 
The argument is based on an observation that in case of tr k < a stable MOTS £ 
with Yamabe type has an intcgrable Jacobi field that leads to a local foliation by 
MOTS on the outside of S, which contradicts the condition that £ be outermost. 
The case of a general tr k can be reduced to this case by bending the data (M, g, k) 
in its ambient space-time to the past, and using the Raychaudhuri equation to show 
that the foliation of MOTS in this new slice gives rise to trapped surfaces outside 
of £ in the original data set. 

Collecting these results, we arrive at the following comprehensive theorem about 
the existence, regularity, and properties of the trapped region [H [31 HH [13 [131 HI] ■ 

Theorem 4.6. Assume that (M,g,k) is an asymptotically flat initial data set of 
dimension 2 < n < 7. There is an explicit constant C > depending only on the 
geometry of (M, g, k) such that the following hold: 

If the trapped region T of (M, g,k) is non-empty, then dT is a smooth, embedded, 
outermost and stable MOTS. The area and the second fundamental form of dT are 
bounded by C and its outward injectivity radius is bounded below by ^ . Furthermore, 
dT is C -almost minimizing with respect to variations in M \ T ■ 

If (M, g, k) satisfies the dominant energy condition, then dT is the union of 
components with non-negative Yamabe-type. If(M,g,k) is a slice of a space-time 
satisfying the dominant energy condition, then the components of dT have positive 
Yamabe-type. 

For the explicit dependence of the constant, see the original references [2l [31 IT51 



To conclude, we wish to point out that the existence of the trapped region in 
(M, <7, k) allows the construction of blow-up solutions to Jang's equation. These are 
nontrivial solutions to Jang's equation which are defined on M\ (7~U f2_) where 
f2_ C M is such that the boundary components of fi_ disjoint from dT arc MITS. 
The construction uses the techniques discussed in section 13.51 and is described in 
[4"3] . A catch however is that some or all of the components of dT may lie in the 
interior of f2_ if they are enclosed by surfaces £ satisfying Hs — trs(fc) = 0. 



In this section we indicate a couple of directions for further research related to 
the ideas discussed in this survey. 

5.1. Generalizations of Jang's equation. The Penrose inequality 



5. Outlook 
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is an equality only for slices in the Schwarzschild spacetime. As we have seen, 
Jang's equation was motivated by the idea of "detecting" data sets which generate 
a Minkowski geometry — dt 2 +g aat . Based on this observation, it appears reasonable 
that any approach to proving the general Penrose inequality must utilize a setting 
which is sensitive to the Schwarzschild geometry. Motivated by this line of thought, 
Bray and Khuri [TTJ[TU] recently extended Jang's equation to a system of equations 
which is designed to identify slices of the Schwarzschild space-time. 

Recall that the Schwarzschild spacetime in isotropic coordinates can be written 
as a warped product with metric g Schw ~ dt 2 (f> 2 where 

The condition that an initial data set (M, g, k) can be represented as a graph 
(x,u(x)) in the Schwarzschild spacetime is then that 

9ab = gah hw ~ (f> 2 D a uD b u, k ab = n ab , 

where Tr ab is the second fundamental form of graph u in the Schwarzschild space- 
time. As shown by Bray and Khuri, one may also in this more general situation 
introduce defects in terms of which the condition that (M,g, k) is the data induced 
on graph(it) in the Schwarzschild spacetime can be characterized. As in the classi- 
cal setup, these data can be calculated in terms of a related Riemannian spacetime, 
which is a product over (M, g). In the generalization this is a warped product over 
M with warping function (j> 2 , i.e. (M x R, g + <f> 2 dt 2 ). This spacetime is additionally 
endowed with a symmetric 2-tensor K which is a lift of k, the second fundamental 
from of M in the spacetime, to the warped product. Recall that in the classical 
Jang equation, the lift of k is simply 7r*fc, where tt is the vertical projection. For 
the generalized Jang's equation, the lift K is defined as 

K = 77* k + (j)d(j}(N)dt 2 

where N is the downward pointing normal of graph(it) in the warped product. 
The generalized Jang's equation now takes the form 

(27) R A ' I -tr Aj K = 0, 

cf. [TU1 section 2]. Due to the lack of symmetry in the warped product, it is 
necessary to consider the warping function <f> as an unknown and add an equation 
for this as well. 

Bray and Khuri [TTJ [TU| propose three different systems of equations incorporat- 
ing the generalized Jang's equation together with equations for 0, which have the 
potential for yielding a proof of a Penrose inequality. As shown by the counter- 
example of Carrasco and Mars [T2] , one version of the Penrose inequality proposed 
by Bray and Khuri, in terms of generalized apparent horizons, is not valid. How- 
ever, in spite of this counter-example, the approach introduced by Bray and Khuri 
may still be applicable to other versions of the Penrose inequality, see the survey 
paper [40] for further discussion. 

The analysis of the systems proposed by Bray and Khuri is made more difficult 
by the fact that <f) tends to zero at the horizon and as a consequence the generalized 
Jang's equation is degenerate there. Bray and Khuri have been able to carry out 
the necessary analysis in the spherically symmetric case, providing a new proof of 
the general PI in this restricted case. 

5.2. Evolution of MOTS. Consider a spacetime which is the maximal develop- 
ment of asymptotically flat data on M for an Einstein-matter system satisfying the 
DEC. Supposing that the Cauchy surface contains a stable MOTS S, which we 
without loss of generality can assume to be outermost, the spacetime contains a 
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black hole, and under some weak genericity conditions the MOTS lies on a spacelike 
marginally outer trapped tube (MOTT) T-L. The MOTT is determined by a choice 
of Cauchy foliation of the spacetime. 

The outermost MOTT is, with the exception of jump-times (see below) space-like 
in the generic case. Thus, the MOTT is an outflow boundary for causal equations 
in its exterior, and the maximal development of the restriction of the Cauchy data 
on M contains the exterior to T~L. 




M 

Figure 3. 

This leads to an exterior Cauchy problem for eg. the Einstein equations in 
spacetime harmonic coordinates. Let a Cauchy surface M be given, containing an 
outermost MOTS. The exterior Cauchy problem is the initial-boundary value prob- 
lem for the evolution of this system in the closed exterior of the MOTT, including 
the MOTS boundary, evolving from the outermost MOTS. This problem can be 
expected to be relevant for the problem of Kerr stability, and in particular it is 
interesting to prove a useful continuation criterion for it. 

If we consider the maximal extension of the MOTT to the future in a spacetime 
with a regular Cauchy foliation, one expects that after a finite sequence of jumps 
[3] this eventually approaches the event horizon. It is an interesting question to 
understand the details of this scenario. In particular, in terms of the Kerr stability 
problem, one expects to have Price law decay of the matter and gravitational energy 
flux across the event horizon. It is reasonable to speculate that the corresponding 
statement holds for the fluxes across the (weakly) spacclike MOTT. As the strength 
of the flux decreases this has the effect of turning the MOTT null. 

This leads to the expectation that the MOTT asymptotically approaches the 
event horizon and terminates at future timelike infinity. Since the MOTT is ex- 
pected to rapidly turn null, one expects the distance along the MOTT to its bound- 
ary at future timelike infinity to be finite. This behavior was verified in the spheri- 
cally symmetric case by Williams (52] who showed that for an Einstein-scalar field 
spacetime with decay along the event horizon of the form v~ 2 ~ s , the MOTT has 
its boundary at a finite distance. As pointed out by Williams, the required decay 
is weaker than the expected Price law decay of v~ 3 . If this scenario is correct, it 
is likely there is a relation between the decay of fluxes across the MOTT and the 
regularity at the boundary of the MOTT at future timclikc infinity. We mention 
here also the work of Ashtekar and Krishnan [SJ [7] in the dynamical horizon (DH) 
setting which shows that the area of the cross sections of a DH is increasing (a 

■*"A MOTT is a dynamical horizon if it is spacelike and foliated by marginally trapped surfaces, 
i.e. MOTS which also have negative expansion with respect to the ingoing null normal, see [5] 
section 2.2] for details. 
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quasi- local version of Hawking's area law for event horizons), and gives expressions 
for the rate of increase of area in terms of the flux across the DH. 

As discussed in section 14731 once a MOTS is created in an evolving spacetime on 
a Cauchy surface Mo then, if the spacetime satisfies the null energy condition, each 
Cauchy slice in the future of Mq contains an outermost MOTS. Further, each time 
a MOTS So is created, it is through a bifurcation process which leads to an inner 
and an outer branch of the MOTT originating at So- The outer branch may jump 
but remains stable, while one expects that the inner branch eventually becomes 
unstable. 

It is of interest to understand in more detail the space-time track of the MOTS. 
The generalized maximum principle for MOTS, cf. section f4. 31 implies that two lo- 
cally outermost MOTS which approach sufficiently closely must eventually be sur- 
rounded by a MOTS. In terms of the evolution of binary black hole data this means 
that two black holes (as determined by their apparent horizons) which approach 
sufficiently closely, eventually are swallowed by a larger black hole surrounding the 
two. 

Ashtckar and Galloway j5j proved a uniqueness result which gives further infor- 
mation on the spacetime geometry of dynamical horizons, a special case of MOTTs. 
This result states that in a spacetime satisfying the null energy condition, the past 
domain of dependence of a DH cannot contain a marginally trapped surface, see [5J 
Theorem 4.1]. It would be interesting to understand better whether results of this 
type hold for MOTTs and MOTSs. 

If one considers two BH's, one of which is small relative to the other, it is natural 
to consider a scenario where the small BH falls into the larger one. In this case, 
the generalized maximum principle for MOTS does not give any information about 
the small BH crossing the horizon of the large one, but the classical maximum 
principle prevents one MOTS from "sliding" inside another. In particular, the 
configuration shown in fig. 15.21 corresponding to the moment when the small BH 
moves inside the larger BH is ruled out by the maximum principle. Therefore 
one expects that as the BH's coalesce, the two apparent horizons will eventually 
approach each other and merge. It is interesting to speculate whether the MOTS 
in such a situation form a continuous spacetime track, with one branch connecting 
the merging horizons with the outermost, surrounding, MOTT. See [44, for details. 

6. Concluding remarks 

In this paper we have given a survey of the state of the art concerning Jang's 
equation, MOTS, implications on the existence of black holes and related issues. 
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The main motivation for considering these issues has so far been in the asymptoti- 
cally flat case. However, it is important to recall that also in considering the Cauchy 
problem for the Einstein equations in strong field situations, analogues of MOTS 
and trapped regions can be expected to play an important role, and therefore some 
of the topics discussed in this survey may have applications in future approaches 
to global evolution problems and the cosmic censorship problem. 
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